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H H ' ratio combmmg (MRC) receivers. The operating environment is Rayleigh-fading with both transmit and 

CJ . receive spatial correlation. We present exact expressions for the probability density function (p.d.f.) of 

the output signal-to-noise ratio (SNR), as well as the system outage probability. The results are based 

, on explicit closed-form expressions which we derive for the p.d.f. and c.d.f. of the maximum eigenvalue 

' of double-correlated complex Wishart matrices. For systems with two antennas at either the transmitter 

, or the receiver, we also derive exact closed-form expressions for the symbol error rate (SER). The new 

' expressions are used to prove that MIMO-MRC achieves the maximum available spatial diversity order, 

c/j , and to demonstrate the effect of spatial correlation. The analysis is validated through comparison with 
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^ ' Monte-Carlo simulations. 
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I. Introduction 

Multiple-input multiple-output (MIMO) antenna technology can provide significant improve- 
ments in capacity [1^] and error performance [5] over conventional single-antenna technology, 
without requiring extra power or bandwidth. When channel knowledge is available at both the 
transmitter and receiver, MIMO transmit beamforming with maximum-ratio combining (MRC) 
receivers [6] is particularly robust against the severe effects of fading. This robustness is achieved 
by steering the transmitted signal along the maximum eigenmode of the MIMO channel, resulting 
in the maximization of the signal-to-noise ratio (SNR) at the MRC output. 

Recently, MIMO-MRC has been investigated in uncorrelated and semi-correlated channel sce- 
narios (i.e. where correlation occurs at only one end of the transmission link, or not at all). A key 
to deriving analytical performance results is to statistically characterize the SNR at the output of 
the MRC combiner. In [7-11], uncorrelated Rayleigh fading was considered, and the output SNR 
statistical properties were derived based on maximum eigenvalue statistics of complex central 
Wishart matrices. In [12], uncorrelated Rician channels were characterized using maximum 
eigenvalue properties of complex noncentral Wishart matrices. Semi-correlated Rayleigh channels 
were considered in [13], utilizing properties of semi-correlated Wishart matrices. 

In this paper we consider double-correlated Rayleigh channels, by first deriving results for 
the eigenvalue statistics of double-correlated complex Wishart matrices. In practice, double- 
correlated channels (i.e. with correlation at both the transmitter and receiver) commonly occur 
due to, for example, insufficient scattering around both the transmit and receive terminals, or to 
closely spaced antennas with respect to the wavelength of the signal. While there are numerous 
statistical results on general Wishart matrices, there are almost no results for the eigenvalue 
statistics in the case of double-correlated Wishart matrices. In [14], the joint probability density 
function (p.d.f.) of the eigenvalues of such matrices was derived in terms of hypergeometric 
functions of three matrix arguments. In [15], the marginal p.d.f. of an arbitrary unordered 
eigenvalue was derived. Here, we add to these general statistical results, by deriving new exact 
closed-form expressions for the p.d.f. and cumulative distribution function (c.d.f.) of the maximum 
eigenvalue of double-correlated complex Wishart matrices. 
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The new general statistical results allow us to consider performance measures for MEMO-MRC 
in double-correlated Rayleigh environments. In particular, we present explicit expressions for the 
p.d.f. of the output SNR, and for the outage probability. These expressions are in closed-form, 
simple, and apply for arbitrary antenna configurations, and over the entire range of SNRs. 

We also derive explicit closed-form expressions for the average symbol error rate (SER) with 
various modulation formats. SER expressions are only currently available in closed-form for 
uncorrelated Rayleigh channels with 2 x m (i.e. 2 transmit and m receive antennas), or m x 2 
systems, where m > 2 [8, 9]; and all other results either require numerical evaluation of unknown 
coefficients (e.g. [7]), or evaluation of infinite series. For the SER results in this paper, we also 
restrict attention to 2 x m and m x 2 systems in order to obtain closed-form solutions, however 
our results are much more general since they apply for double-correlated channels. The 2- 
antenna restriction in fact has important practical applications. As discussed in [9], size and cost 
constraints in realistic cellular MIMO systems will typically limit mobile units to accommodate 
a maximum of two antennas. Base stations, on the other hand, usually have greater design 
flexibility, and can easily be implemented with more antennas. As such, 2 x m and m x 2 
correspond to uplink and downlink communication scenarios in such systems, respectively. 

Even more insights can be gained by analyzing the performance at high SNR. In particular, 
we formally prove that MIMO-MRC achieves the maximum possible spatial diversity order in 
double-correlated Rayleigh environments, and show a clear and direct relationship between the 
SER and the transmit and receive correlation matrices in the high SNR regime. 

Finally, we verify the analytical SNR p.d.f., outage probability, and SER expressions by 
comparing with Monte-Carlo simulations, and examine the impact of correlation in each case. 
Our results show that the outage probability may increase or decrease as a function of the spatial 
correlation, whereas the SER tends to vary monotonically. 

II. Statistical Properties of Double-Correlated Complex Wishart Matrices 

This section presents new general statistical properties of complex double-correlated Wishart 
matrices, which will be used in the subsequent performance analysis of MIMO-MRC. 



Submitted to the IEEE Transactions on Communications (Submitted November 2005) 



3 



A. Cumulative Distribution Function ( C.D.F. ) 

The following theorem presents the c.d.f. of the maximum eigenvalue of double-correlated 
complex Wishart matrices. This will be used for deriving the outage probability of MIMO-MRC 
in double-correlated Rayleigh channels. 

Theorem 1: Let X ~ CMm,n (Omxn, S ® U), where n < m, and O G C'^^'^ and S e C'"^"* 
are Hermitian positive-definite matrices with eigenvalues Ui < . . . < Un and Ui < . . . < Um 
respectively. Then the c.d.f. of the maximum eigenvalue of the double-correlated complex 
Wishart matrix X^X is given by 

(-l)"r„Hdet(f^)"-Met(Er-Met(^(a:)) 

A„(rj)A„(S)(-a;)-("-i)/2 ^ ^ 

where r„(-) is the normalized complex multivariate gamma function, defined as^ 

n 

r„(n) = J]r(n-z + l) (2) 

i=l 

and A^( ) is a Vandermonde determinant in the eigenvalues of the m-dimensional matrix 
argument^, given by 

m 

A^(S) = J]((7,- - a,) . (3) 

i<j 

Also, '*^{x) is an m X m matrix with element 

{/ \ m—i 
^" 'L . . (4) 

e "i-r-ip 1^^; for I > r 

where t = m — n, and 

V{t,y) = l-e-yY.y^/k\ (5) 

fe=0 

is the regularized lower incomplete gamma function. 

Proof: See the Appendix. ■ 

'Note that this is related to the standard complex multivariate gamma function r„(n) (as defined in [16]) via r„(n) = 
7i--"("-i)/2f„(n). 

^For notational convenience, we will sometimes give a set in place of a matrix in the argument of Am(-), as a shorthand. 
Specifically, for the set A, we use A^(^) as shorthand for A,„(diag(.A)). 
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Corollary 1: For the case n = 2,m > 2, ^ reduces to 



det (n) " 



^^-^"^^ = A (17) A (S) 5Z 5Z (-1)"^'^*^ M"^-" (al^'*]) Q,,(x) (6) 



where 



t ,t<p 

0W = <! (7) 

t - 1 , t > p 
and = {a,; i e {1, . . . ,m}\{p, t}}, and 

1 ^„ / X \ 3^„ / X 



Qptix) = -e rn; e "i-tP m; . (8) 

Proof: The proof follows by applying Laplace's Expansion to the last two rows of det {^{x) ) 
in ([T]), and then simplifying the resulting (Vandermonde-type) minors using (15^ . ■ 
Corollary 2: For the case n = m = 2, ^ reduces to the simple expression 



Proof: The proof is straightforward and is omitted. 



B. Probability Density Function (P.D.F.) 

The following theorem presents the p.d.f. of the maximum eigenvalue of double-correlated 
complex Wishart matrices. This will be used for deriving the p.d.f. of the output SNR of MIMO- 
MRC in double-correlated Rayleigh channels. 

Theorem 2: Let X ~ C7V„,„ (0„xn, S ® ft), where n < m, and CI e C"^" and S G C™^"^ 
are Hermitian positive-definite matrices with eigenvalues cui < . . . < cUn and ai < . . . < am 
respectively. Then the p.d.f. of the maximum eigenvalue of the double-correlated complex 
Wishart matrix X^^X is given by 



.(.-l)detM ^^^^^^^^^^^ ^ ^^^^ 



2A.,. 

where ^^(Am) is an m x m matrix with element 



= T+1 
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:*.(A.)),, = \ . (11) 

" ' vim- 1; — ^ ) foTi = i 



and where {^{Xm))^^ is defined in ©. 

Proof: The result follows by differentiating (HI) with respect to x. ■ 

III. MIMO-MRC System Model 

Consider A^^ transmit and A^,. receive antennas, where the A^^ x 1 received vector is 

r = Hwx + n (12) 

where x is the transmitted symbol with = 1, w is the beamforming vector (specified 

below) with E [||w||^] = 1, n is noise ~ CM Nr,i (Oat^xi, '^Nr)^ 7 is the transmit SNR. Also, 
H is the Nr x A^^ channel matrix, assumed to be flat spatially-correlated Rayleigh fading, and is 
decomposed according to the common kronecker structure (as in [4, 14, 17, 18], among others) 
as 

H = R^H^S5 ~ CAr^^,^,(0^,x^„R®S) (13) 

where R and S are the receive and transmit correlation matrices respectively, with unit diagonal 
entries, and ~ CM NrM (Ojv.xiVt, I^v. ® I^vJ- 
The receiver employs the principle of MRC to give 



z 



w"^HV = Vfw'^H'^Hwx + w"fH"^n . (14) 



Therefore, the SNR at the output of the combiner is easily derived as 

7 = 7wWHw . (15) 

The BF vector w is chosen to maximize this instantaneous output SNR, thereby minimizing 
the error probability. It is well known that the optimum BF vector Wopt is the eigenvector 
corresponding to the maximum eigenvalue \m of H^^H. In this case, the output SNR dTSt becomes 

7 = 7wlptH"^Hwopt = l\m ■ (16) 
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Clearly the output SNR (and therefore the performance) of MIMO-MRC depends directly on 
the statistical properties of Am- 

In the previous section we presented results which apply directly to Am for the case when 
Nr > Nt (i.e. H^H is full-rank), by setting fl = S, T, = R, n = Nt, and m = Nr. The resuhs 
of the previous section also apply to the case when Nr < Nt, (i.e. HH^^ is full-rank), by setting 
17 = R, S = S, n = Nr, and m = A^^; since the maximum eigenvalue is the same in both cases. 

IV. Performance Analysis of MIMO-MRC in Double-Correlated Channels 

A. Statistical Characterization of the Output SNR 

We now characterize the statistics of the output SNR 7. Using (O, and making a simple 
change of variables to (flOb we obtain the p.d.f. of 7, given in closed-form by 



The outage probability is an important quality of service measure, defined as the probability 
that 7 drops below an acceptable SNR threshold 7th. It is obtained from ([T]) and (fTSl) as follows 



In Section |Vl we show that the outage probability can increase or decrease with the spatial 
correlation, depending on the value of 7th- 

B. Symbol Error Rate Analysis 

In this subsection we derive closed-form expressions for the average symbol error rate (SER) 
of MIMO-MRC in double-correlated Rayleigh channels, with various modulation formats. As 
discussed in Section ID we consider the practical special cases 2 x m and m x 2. Our results 




(17) 




(18) 
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apply for all general modulation formats that have a SER expression of the form 



(19) 



where Q(-) is the Gaussian Q-function, and a and b are modulation- specific constants. Such 
modulation formats include BPSK (a = 1, 6 = 1); BFSK with orthogonal signalling (a = 1, 6 = 
0.5) or minimum correlation (a = 1, 6 = 0.715); and M-ary PAM (a = 2(M - 1)/M, b = 
3/(M^ — I)). Our results also provide the approximate SER for those other formats for which 
ST% is an approximation, e.g. A//-ary PSK {a = 2,b = sin^(7r/M)) [19, Eq. 5.2-61]. 

A common useful approach for evaluating SERs of the form (fT^ . is to first evaluate the moment 
generating function (m.g.f.) of 7, and then apply the well-known m.g.f.-SER relationships given 
in [20]. In the context of MIMO-MRC, this approach was used in [7,9] to evaluate SERs in 
uncorrelated Rayleigh channels. Although it is possible to evaluate the output SNR m.g.f. for 
double-correlated channels considered in this paper, unfortunately the resulting expression has 
significant convergence problems. 

As such, we adopt an alternative approach. We begin by generalizing a result presented in 
[13], which directly related the SER of MIMO-MRC in semi-correlated channels with BPSK to 
the c.d.f. of the output SNR. We start by noting that the Gaussian Q-function can be written as 



Qix) 



/V2 



and substitute into (fT9t to give 

P. = 



dv 



bu 



f-y{u)du . 



(20) 



(21) 



Now we apply (definite) integration by parts to d^TTl . To do this we integrate fj{u), and differ- 
entiate the quantity in square brackets as follows 



d_ 

du 



e-'" dv 



bu 



d_ 

du 
d_ 
du 

1 lb 



Vbu 



e"" dv 



d_ 

dx 



e-" dv 



x=\/bu 



-bu 



(22) 
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where the first line followed by noting that 



n 



(23) 



and the second line followed from the Chain Rule, and last line was obtained using the Second 
Fundamental Theorem of Calculus. This gives 

a\/b 



P. 



oo g-felt 



^ , ^-Fy (u) du. 



(24) 



Note that (l24ll is a generalization of [13, Eq. 32] to SER expressions of the form (fT9l) . 

Now we can directly substitute (fTSt into (l24ll and integrate term by term. This procedure, 
however, is problematic since the term-by-term integrals are divergent (although the overall 
result converges). To circumvent this problem, we first write (l24l) in terms of the c.d.f. of Am, 
given in Corollary [H as follows 



oo g-fell 



7 



2aA Jo a/u 
and derive an alternative expression for Qp^t{x) in dSj). 
Expanding the exponentials into power series, and using we manipulate dH) as follows 



(25) 



1 

X 



k=m 



UJ2<Jp 



k\ 



( 



k=m 



\ 



k\ 



oo oo ( r 



\ki+k2 j 1 



ci=m 
oo 

E 



1 

UJlUt 



k2 



ki=m k2=m 

oo 

/ \k—l (-> 



(26) 



where 



k=2m 



i=m ^ ' 



k-e 



Now, substituting (l26b into Q and (l25t . the SER can be written as 

aVb det (O) 



m m 



20FA2 (fl) A, 



E E( 

p=l t=l,i7^p 



(27) 



(28) 
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where 



k=2m 

CO / 

= 7V^5^5J--J T{k-l/2) (29) 

k=2m ^ '^^ 

where the last line followed by applying the integration identity [21, Eq. 3.381.4]. To remove 
the infinite summation from d^^ . we apply the Binomial Theorem to express Sk in (ITTt as 
1 // 1 1 \' 'izi/k\ / 1 \'/ 1 



Now, performing the change of variables i {k — t) va the last summation in (l30b . and using 



the property (^) = we can write 



1(1 1 \ s-^ I k\ ( 1 \ ( 1 \ / 1 \ / 1 



(31) 

Note that the key advantage of writing Sk in this way, as opposed to (l27t . is that we have 
removed the dependence of the summation limits on k. We now substitute dSTl) into ( l29t and 
simplify to obtain 



1 \\ V ■\ 1 ^ ( ' 760^2 CTp ) V ■\ 1 l^' ■ybuJiat 



i 

pj 



(32) 



where 



fc=2m 

r (fc + £ - 1/2) 



E 



A;! 

fc=2m-f 



r a - 1/2) ,Fo (£ - 1/2; y) - '^J^' ^ + f '^'^ (33) 



fe=0 
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where iFq (■) is the binomial hypergeometric function. Noting that iFq (a; y) = 2F1 (a; 6; 6; y), 
and using [22, Eq. 15.1.8], we obtain 



k\ 

fc=0 



(34) 



Now we apply the identity 



where 



r(,.l/2).2t^ (35) 



(2A;- 1)!! = 1 X 3 X ... X (2fc- 1) (36) 



for /c > 0, and define (—1)!! = 1 and (—3)!! = —1, which yields 



2« 

where 



v[i'.y) = ^,W,v) (37) 



k=0 

Finally, substituting (1771) into and ( ESI) yields the desired closed-form SER expression 

' ™ ^ ^ p=i t=i,i^p 

Note that (l39b is a simple finite closed-form expression (involving only polynomial terms in 
7), which can be evaluated easily and efficiently. 

C. High SNR SER Analysis 

We now analyze the SER performance in the high SNR regime in order to derive the diversity 
order of the system. 

Consider the SER expression given by (l28t and (l29t . It can easily be shown that, irrespective of 
the SNR, the first term in (i.e. k = 2m) cancels with terms in ( ESI to give zero contribution 
to Ps- Therefore the summation in (j29|) could equally well be written starting from k = 2m + 1. 
Since as 7 — 00, the infinite series in (E^ is dominated by the low order terms, in the high 
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SNR regime we need only to consider k = 2m + 1. To proceed, we take the k = 2m + 1 term 

in the series and use (l27t to write ( l29t at high SNR as follows 

^ (7&)-^-r(2m + l/2) ((J_Y (J_Y^' ^ /J_Y^' 

\/bm\{m + l)\ \\u;20-pj \uJiat J \uJ2crpJ V^iO"t. 

(76)-'™ r(2m + l/2) 1 f 1 1 



Vbm\{m + 1)1 det (f2)™ (o-p(7t)™ V^i*^* ^2(7^, ' ^^^'^ 



Substituting dJOl) into (EH) gives Ps at high SNR as follows 

ar(2m+ 1/2) (76)""" 
2v/?m!(m + 1)!A2 (O) A„ (S) det {ft 

where 



poo ^ -^y-'- ' ^ (41^ 



m m 



S = y y (^ + ^] ■ (42) 



Now, it can be shown that 



m m 



^ ^ and ^ (-ir^W ^-\^^'"'^^ 

are Laplace expansions of (— l)™^^A„i(S) and (— l)™~^Am(S) respectively. As such, (H^J can 
be written as 

^ ,1 1 \ A^fSl 



c^i 7 det (S)^ 
^ A2(^2)A^(S) 

det (fl) det (S)^ 

where the second line followed from (l56t . Finally, substituting (l43t into (HTl) and simplifying 

using (l35l . we obtain the desired high SNR SER result 

poo ^ a (4m -1)!! . 

52m 22m+i ^!(^ + 1)1 det (ri)" det (S)^ 
Therefore, clearly MIMO-MRC achieves the maximum possible spatial diversity order of 2m in 

double-correlated Rayleigh channels. Also, using Hadamard's inequality [23, Th. 16.8.2], and 

the fact that the diagonal elements of Ct and S are unity, it can be easily shown that 

< det (O) < 1 and < det (E) < 1 (45) 

with equality in the upper limit only when the correlation matrices are identity matrices. This 
proves that the presence of spatial correlation (at either end) yields a net reduction in error 
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performance in the high SNR regime. 

V. Numerical Results 

While the analytic results in this paper apply to arbitrary channel correlation matrices, for 
our numerical studies we construct the correlation matrices using the practical channel model 
presented in [24]. The model assumes that there are uniform linear arrays at both the transmitter 
and receiver. Let us denote the relative antenna spacing between adjacent antennas (measured in 
number of wavelengths) as d,. at the receiver and dt at the transmitter. Also, define 9r, 9t, and 
cr^ as the mean angle of arrival (AoA), mean angle of departure (AoD), receive angle spread 
and transmit angle spread respectively, and let 6_j. = 6^. + Or and 9t = Ot + dt denote the actual 
AoA and AoD, with 6,. ~ M (0, a^) and Ot ^ M (0, a"^). With these definitions, the (p, qf"^ entry 
of R and S is given by 

R_ -j2Tx{q-p)dr cosiOr) ^-h{2iv(q-p)dr s\n{er)crr)'^ Q _ -j27r(p-i})dt cos(6't) -i(27r(p-(j)di sin{et)o-i)^ 
p,q — a a ^ J 'Jp,q — a a ^ 

Recall that R and S are directly related to O and E as discussed at the end of Section [nil For 
all results in this section we assume = dt = \ and 6^- = 6t = ^■ 

Fig. [T] shows the p.d.f. of the output SNR of MIMO-MRC with various antenna configurations. 
The 'Analytical' curves are from (flTt . and clearly agree with the Monte-Carlo simulated p.d.f.s. 
Moreover, we observe that both the mean and variance of the output SNR increase with the 
number of antennas. 

Fig. 12 shows the (analytical) p.d.f. of the output SNR, comparing various correlation scenarios. 
We see that for 4x4 antennas, extra correlation increases the spread of the SNR around the 
mean. This effect was also observed for semi-correlated channels in [13]. For the 2x2 case, 
the correlation has less effect. 

Fig. |3] shows the outage probability of MIMO-MRC, with the same antenna configurations as 
in Fig.[TJ The 'Analytical' curves are from dTSt . and agree precisely with Monte-Carlo simulated 
curves. We see that the outage probability is significantly improved as the number of antennas 
are increased. 

Fig. m shows (analytical) outage probability curves, comparing different correlation scenarios. 
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We see that for both antenna configurations, the correlation increases the outage probability at low 
SNR thresholds, and decreases the outage probability (thereby improving system performance) 
at high SNR thresholds. These general trends were also previously observed for semi-correlated 
channels in [13]. Moreover, we see that the cross-over point of the different correlated curves 
occurs at lower outage probabilities as the numbers of antennas increase. 

Fig. 13 shows the SER of MIMO-MRC with BPSK modulation, for various antenna configu- 
rations. The 'Analytical' curves are from with a = 1 and b = 1, and match exactly with 
the Monte-Carlo simulated SERs. The 'Analytical (High SNR)' curves are from (H?b . Clearly 
they converge to the exact SER in the high SNR regime, confirming that the maximum possible 
diversity order is achieved. 

Fig. |6l shows the SER of MIMO-MRC with 4-R\M modulation, for various antenna configu- 
rations. The 'Analytical' curves were generated from (l39l) with a = 1.5 and b = 0.2. Again we 
see an exact agreement with the Monte-Carlo simulated curves in all cases. 

Fig. |7] shows the SER of MIMO-MRC with QPSK modulation, for various antenna config- 
urations. The 'Analytical' curves were generated based on (l39b with a = 2 and b = 0.5. As 
discussed in Section IIV-BI (l39b only provides an approximation for QPSK, however we see that 
this approximation is accurate for all but very low SNRs. In particular, for all SERs of practical 
interest (i.e. P, < 0.01), the analytical curves match almost exactly with the simulated curves. 

Fig. |5] uses the analytical expressions to examine the effect of correlation on the SER. Results 
are presented for BPSK modulation. We see that for both antenna configurations, the SER 
increases monotonically with the level of correlation. This agrees with the high SNR predictions 
from dUl) and ( B31 . 

VI. Conclusion 

We have examined the performance of MIMO-MRC in double-correlated Rayleigh fading 
environments. Our results are based on exact closed-form expressions which we derived for the 
p.d.f. and c.d.f. of the maximum eigenvalue of double-correlated complex Wishart matrices. We 
showed that the outage performance may increase or decrease due to the presence of spatial 
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correlation, depending on the average SNR. When the minimum number of transmit and receive 
antennas is two, we proved that MIMO-MRC achieves the maximum available spatial diversity 
order in double-correlated channels. 

Appendix 

A. Proof of Theorem El 

First consider the n = m case of square random matrices X ~ CM m,m{^mxm, S ® f2), and 
let Ai < . . . < Am be the non-zero eigenvalues of X^^X. The c.d.f. of A^ is obtained using 

Fx^{^)= f /(A)rfA (46) 



where A = diag{Ai, . . . , Am}, /(A) is the joint p.d.f. of Ai, . . . , Am, and I? = {0 < Ai < . . . < A^ < x}. 
It was shown in [14] that 

^ oFo(-fl-M]-\A)Am(A)^ 
■'^ ' rm(m)2 det(^^)-det(S)- ^ ^ 

where o-Po (■;■;■) is a complex hypergeometric function of three matrix arguments. To evaluate 
the integral in we first expand o-fo(') in complex zonal polynomials [16] as 

„F„(-n-\S-\A)^f y:^^<-"''),^^<^''fe<^' (48) 

fc=0 K. fi-.0/cl,lmJ 

where the inner sum is over all partitions /C = (/ci, . . . , km) with /ci > /i;2 > • • • > 0, and 
fci + . . . + fcm = ^- Using the character representation for complex zonal polynomials [16] 

C^(A) = X[^](1)XW(A) (49) 

and 

~ X[/c](l)'rm(m,/C) 

where is the character of the representation {/C} of the linear group, X[/c](l) is the 

dimension of the representation [/C] of the symmetric group, and 

rm(m,/C) = JJr(m + A;i-^+l) (51) 
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we can write 



oF„ (-n-.s-. A) , 5^4^)^^q(-n-')x^x,(I:-')x^x,(A) ^j^) 



fc=0 IC 

Now we apply Weyl's formulas [25] 



'I) = k\ 



(53) 



and^ 



det a: 



X{x;}(A) 



kj+m-j 



det (A--^) 



det a: 



kj+m-j 



A„ (A) 



(54) 



to give 



oFo(-ri-\S-\A) 



n 



i<j 



det 



EE 



------ 



det ^ 



kj+m—j 



det a: 



kj+m-j 



(55) 



fc=0 /c 



rK/c)n™,(fci-fc.-j + OA„(A) 

Now changing from /C to strictly ordered partitions ICo = (ki, . . . ,km) with fci > . . . > fc^ > 
(i.e. such that ki + m — i ^ ki), and using 

1 1 \ _ UT<j i^i - ^j) 



n 

Kj 



UJj UJi 



we obtain 



nm m-1 
i=l 

m—1 1 , /■^\m—l 



(56) 



(-l)"(™-i)/2r^(m)2 det (17)'"^' det (S 



A™(rj) A„(5]) 



X 



EE 



det 



det 



det A 



Substituting (l57t into (l47t and simplifying, the joint eigenvalue p.d.f. becomes 

_]^^m(m-l)/2 



(57) 



/(A) 



det(rj) det(S)A^(ri)A^(S) 



EE 

A;=0 /Co 



det 



det 



det A^ A™(A) 



(58) 



^Here we introduce the compact notation for the determinant of a matrix, written in terms of the {i,jy^ element. 
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Now substituting dSSt into (l46b . the c.d.f. of can be written as 



where 



-1) 



m,{m—l)/2 



det 



det(rj) det(S)A„(ri)A„(S; 



EE 

fc=0 ICo 



det 



J 



(59) 



J= / det AfM A„(A)rfA. 



(60) 



To evaluate this integral, we note that (A) = det (A] ), and apply [26, Corr. 2] to obtain 



J = det 



m(m+l) 
X 2 



det 



(61) 



Now, in order to remove the infinite sum over partitions in (l59l) . we must manipulate the right- 
hand determinant in (|6T1) as follows 



det 



1 



kj + i 



det 



fcl+l fc2+l 

1 1 
fcl+2 fc2+2 



'^77^ + l 

1 

fcm + 2 



\ fcl4 



1 



(62) 



det 



fcl— fc2 



fcl — fcrf 



fcl+l 
1 


(fcl+l)(fc2+l) 
fcl — fc2 


(fcl+l)(fc„ 
fcl — fcrt 


,+1) 


fcl+2 


(fcl+2)(fc2+2) 


(fcl+2)(fe„ 


^+2) 


1 


fcl— fc2 


fcl — fc. 




fcl +m 


(fcl4-m)(fc2+m) 


(fcl +m) (fc77 





(63) 



11(^1 - ^ ) n 



i=2 



1 \ki+i 



det 



11(^1 - ^ ) n 



i=2 



\ \ki+i 



det 



/l ^ 

fe2 + l 



1 

/ 1 






fcm +1 
1 

fem+2 



(64) 



n+m / 



1 

fc2 + l 
~1 



fcm + 1 
-1 



(fc2+l)(fc2+2) 
-2 

(fc2 + l)(fe2+3) 



(fcm + l)(fem+2) 

^2 

(fc„ + l)(fe„+3) 



(fe2 + l)(fc2+™) 



-("I-l) 



(65) 



(fcm+l)(fem+m) 
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m 

r(m)n 

i=2 



ki ki 
ki + 



+ i 



det 



/ ^ ... 

k2+2 

1 

k2+3 



\ 



fem+3 



(66) 



Note that in (l63t and (l64b we have subtracted the first column from all other columns and removed 
common factors. In (l65l) and (l66b we have then subtracted the first row in the determinant from 
all other rows and removed common factors. Notice that the determinant in (l66t is a principle 
submatrix of the determinant in (1^^ . Applying the same sequence of operations (m — 1 times) 
to the determinant in (l66b . we obtain 



det 



1 \ki + i 



Substituting (l67t into (l6TT) and (l59b and simplifying gives 
ij L^[m)x 5^ 5^ det 



where 



det(f^) det(S) Am{n) A^(S) ^ ^ 



1 



det 



1 



Now we apply the Cauchy-Binet formula [27] to give 

(-l)'"("^-i)/2r„(m)a;™("'+^)/2 
^"'"^"'^ " det(f^)det(S)A^(f^) A^(S) 

and deal with the infinite sum as follows 

k 



det K] 



vfc=0 



CO / -, \ fe oo 
fc=0 ^ '^^'^^^ fc=0 



fc=0 

oo 

E 

fc=0 

1 

m! 



A;! 



r(A; + l) 



k\ 



iFi l;m + l; 



r(m + A; + l^ 

a; 



(67) 



i=l 



(68) 



(69) 



(70) 



(71) 
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where i-Fi(-) is the confluent hypergeometric function. Applying the identity (see [22, Eq. 6.5.2] 
and [22, Eq. 6.5.12]) 

iFi (1; b; z) = V{h)z'-'e' V {b - 1; z) (72) 
in dTTT) . and then substituting into (ITOb yields 

This establishes the result for square matrices. 

To obtain the result for rectangular matrices (i.e. for n < m, and n x n matrix O), we follow 
an approach of [15], and consider an auxiliary (expanded) system with m x m matrix CI with 
eigenvalues (cDi, . . . , ujm) = (ei, . . . , e,-, wi, . . . , u;„), for which dTSt holds, and take limits of (1731 
as ei — > 0, . . . , e^- — > 0. We then have 



where 



_ \ m— 1 

det (O 



To take these limits, we start by noting that 



lim ^ /_^- det ( e "^"iP ( m; --r^ ) ) . (75) 



m; = e "''"j 1 ^ — -. — — ^ 

V y r(m) 



(76) 



where r(-, ■) is the upper incomplete gamma function, and use [22, Eq. 6.5.32] to obtain 

Coiaj) r(m) V Cbiaj 



— _ , X \ — 1/ X ^ f m — 1 (m — l)(m — 2) 
hm e P m: = — - 1 + — + ^ . ^1 — ^ + . . . 



X \2 



(77) 

We then consider each of the r limits in turn, starting with uji ^ 0. We replace the first 
row in the right-hand determinant in (1751) with the first term in the expansion in dTTl (i.e. 

(\ m— 1 
— ) ). For the denominator we have 

m 

lim Am(Sl) = A„_i(u}2, • • • FT '^i • C^S) 
Now, factoring the numerator determinant yields a finite ratio for C in (1751 in the limit cJi ^ 0, 
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given by 



(\ m— 2 
U'^=2^j) det(Hi(x)) 



C = ' '. , lim ^ ^ (79) 

Tim) (1.2^0,. ...lir^o A„_i(u;2, . . . , ujm) 



where S^(x) is an m x m matrix with element 



, m—i 

1 



i^))^.= { ^I'i / X (80) 



e ^i^iV m: A> 



UJiO-j 

Now we take the limit as — > 0. In this case the leading divergence in the numerator comes 
from the second term in the expansion (l77t (i.e. taking the first term results in a determinant 

I 1)771-1 / \m-2 

of zero). Hence, we can replace the second row of Si(a;) with r(m-i) ia^) ■ 
denominator we have 

m 

lim A„_i(0) = Am-2{i^3, ■ ■ .,oJm)Y\'^j ■ (81) 

j=3 

Hence, factoring the numerator determinant again yields a finite ratio for £ as cJ2 ^ 0, given by 

/I \m+(m-l) (m,-l)+(m-2) ( '^i ) det {'B,2{x)) 

^ = ^ lim ^-T 4 ^ . (82) 

T[m)T[m — 1) lj3^o,...,iI>t^o Am-2(i^3j ■ ■ ■ 



Continuing this procedure until all limits have been calculated, and simplifying the result, yields 

^ (-l)^x-^det(ll)-Met(5.(x)) 

r(m)...r(m-r + l) A„(f2) ' ^ ^ 

Now noting that Ht-(x) = ^'(a;), where ^{x) is defined in the theorem, and substituting ( 15^ 
into (17^ . we perform some basic manipulations to obtain the desired result in (d). □ 
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Fig. 1. P.d.f. of the output SNR of MIMO-MRC in double-correlated Rayleigh channels, for 7 = dB. Correlation parameters 
are 61. = St = f , d = i, = ^, and a? = ^. 
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Fig. 2. 



P.d.f. of the output SNR of MIMO-MRC in various double-correlated Rayleigh channels, for 7 = dB. 
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Fig. 4. Outage probability of MIMO-MRC in various double-correlated Rayleigh channels, for 7 = dB. 
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Fig. 5. Symbol Error Rate of MIMO-MRC in double-correlated Rayleigh channels with BPSK modulation. Correlation 
parameters sae Or = 6t = ^, d = 5, = f^, and cr? = 55. 
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Fig. 6. Symbol Error Rate of MIMO-MRC in double-correlated Rayleigh channels with 4-R\M modulation. Correlation 
parameters are = = f , d = |, cr^ = j^, and erf = ^. 
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Fig. 7. Symbol Error Rate of MIMO-MRC in double-correlated Rayleigh channels with QPSK modulation. Correlation 
parameters eae Or = 6t = ^, d = i, cr^ = j^, and cr? = ^. 
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Fig. 8. Symbol Error Rate of MIMO-MRC in various double-correlated Rayleigh channels with BPSK modulation. 



